Ultracold dipolar few-boson ensembles in a triple well trap 

Budhaditya Chatterjee, 1 '^ Ioannis Brouzos, 1 ^ Lushuai Cao, 1 ' 2 '^ and Peter Schmelcher 1,2 '^ 

1 Zentrum fur Optische Quantentechnologien, Universitat Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 

2 The Hamburg Centre for Ultrafast Imaging, 
Luruper Chaussee 149, 22761 Hamburg, Germany 
(Dated: December 20, 2012) 

Abstract 

We investigate the ground state properties and tunneling dynamics of ultracold dipolar bosons 
in a one dimensional triple well trap from a few-body ab-initio perspective. Our focus is primarily 
on the distinctive features of dipolar bosons compared to the contact interacting bosons. Forma- 
tion of intra-well localization is observed for very strong dipolar interaction. General population 
rearangement as well as fragmentation and localization effects have been found, depending strongly 
on the particle number. The energy spectrum for two particles exhibits avoided crossings that 
lead to several distinct resonances involving different bands, i.e. to an inter-band resonant tun- 
neling dynamics. The corresponding mechanisms are investigated by studying among others the 
pair-probability and performing an eigenstate analysis. 
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I. INTRODUCTION 



Recently dipolar quantum gasesJTJ [2] have become a major focus in the field of ultracold 
atoms. The important experimental progress of achieving Bose-Einstein condensates of 
chromium atoms [HI H] and more recently of dysprosium atoms [5], as well as a high phase- 
space density gas of polar KRb molecule [6J has strongly motivated these studies. The 
important feature of dipolar gases is their long range and anisotropic interaction which 
results in interesting effects such as an elongation of the condensate along the direction of the 
oriented dipoles [7HH] or the stability of pure dipolar condensates in special geometries such 
as a very oblate trap [7\ EJ ITUHT2"] . Dipolar gases thus represent an important tool to study 
and simulate a variety of novel quantum effects and exotic quantum phases (2J IT3HT5] .Qf 
special importance is the study of dipolar gases in lower dimensions which often results in new 
phases and effects not observed in three dimensions such as p-wave Fermi superfluids in two 
dimensions [16, 17J and Luttinger-liquid-like behavior in one dimension [T8H21] . Moreover 
dipoles, polarized perpendicularly in one-dimension, interact with repulsive interaction thus 
alleviating the collisional instability arising from the head-to-tail alignment of the dipoles. 

Investigations of dipolar gases in one dimensional lattices demonstrate interesting features 
not observed for bosons with on-site interaction, such as density waves [TBI [22] and Haldane 
insulating phases |22| 123] . Moreover, dipolar gases in a ring geometry have been shown to 
exhibit anisotropic interactions leading to the formation of a variety of new phases [24, 25J. 
Dipolar bosons in a triple well potential have been studied in ref |2"T] showing the non- 
local character of the interaction as well as dynamical superpositions [26J. Also in the case of 
the corresponding quantum dynamics, dipolar gases show a rich behavior as a consequence 
of their long-range nature. For instance, the tunneling dynamics of dipolar degenerate 
gases in double-well trap exhibits symmetry breaking and self-trapping [28J while interesting 
tunneling mechanisms have been reported for toroidal traps [29J . 

In this work, we investigate the ground state properties and tunneling dynamics of ul- 
tracold few dipolar bosons in a one dimensional triple well trap from a bottom up ab- 
initio perspective using the numerically exact Multi-Configuration Time Dependent Hartree 
(MCTDH) method [30, 3 lj . The ground state properties along the crossover from weak to 
strong interactions are studied with a focus on microscopic quantum effects and mechanisms 
occurring in the system. We elaborate the differences of dipolarly versus contact interacting 
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bosonic systems [32]. Fragmentation and localization scenarios for very strong interactions 
as well as inter-well population realignment depending strongly on the filling fractions are 
derived. The energy spectrum for two particles show avoided crossings which allow for in- 
teresting inter-band tunneling processes reminiscent of those for contact bosons in a triple 
well [33]. 

The paper is organized as follows. In Section [IT] we introduce our model and the system 



set up. In Sec. Ill we investigate the ground state for different particle numbers. The energy 



spectrum for two bosons is investigated in Sec. IV and subsequently in Sec. M we explore 



the inter-band tunneling dynamics for this case. Sec. VI contains the conclusions. 



II. MODEL AND SETUP 



Our system consists of N interacting dipolar bosons in a ID triple well trap. The general 
Hamiltonian can be written as 

TV *2 N i 

H = E ~zM%i + ? + 2 E v U*i - (!) 

i=l i=l i^j 

Here, the V tw represents the triple well potential which we model as V tw = Vosin 2 (Kx) 
(where Vq gives the lattice depth and n is the wave vector). In order to confine the bosons 
in the triple well we impose a hard wall boundary condition at x = ±37r/2K. 

The interaction potential Vi n t consists in general of a combination of short ranged contact 
interaction and the dipolar interaction. In order to realize pure dipolar effects, we omit the 
short-ranged contact interactions (they can be experimentally tuned to zero using magnetic 
or optical Feshbach resonances [31]). The repulsive dipolar interaction between two bosons 
whose dipole moments are aligned by an electric field is given by Vd(xi — Xj) = D 2 /\xi — Xj\ 3 
where D 2 = d 2 m j^€Q for electric dipoles and D 2 = d^Ho/^ for magnetic dipoles (d m being 
the dipole moment). 

For reasons of universality as well as computational convenience, we will rescale Eq. [T]in 
units of the recoil energy Er = h 2 n 2 /2M thereby setting % = M = k = 1. The scaled triple 
well potential is given as V^ w = Vosin 2 x with hard wall boundary conditions at x = ±37r/2. 
The scaled dipolar interaction strength becomes D 2 /Ti 2 k = d. The length and time units 
are given in terms of and H/Er, respectively, d thus represents a dimensionless dipolar 
interaction strength. All quantities are in dimensionless units throughout. 
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III. GROUND STATE PROPERTIES 



A. Commensurate Filling 

The ground state configuration is determined by the competition between the external 
trapping confinement and the interaction between the particles. In the case of commen- 
surate filling, the most important effect is the superfluid to Mott-insulator transition. For 
small interactions, the effects of dipolar interactions are almost identical to that of contact 
interactions and it is only for significantly stronger interaction strengths that the two could 
be distinguished, in particular concerning on-site effects. 

As we observe in Fig. [I] for six atoms (filling factor v — 2) increasing the dipolar inter- 
actions leads to an equal population distribution among all the wells. This is due to the 
tendency to minimize the energy since multiple occupation of a single well is energetically 
unfavorable for stronger interactions. The system thus forms a Mott-Insulator with each 
well occupied by two particles. On the top of that, since each well has a double occupancy, 
the density profile in each site is strongly influenced by the repulsive dipolar forces and 
the atoms tend to avoid each other on the same site, indicating a localization of the two 
bosons in each well forming a "crystal-like state". This effect is even stronger for dipolar 
interactions than for contact repulsion [32j due to their long-range nature. While in the case 
of contact interacting bosons, the interaction energy saturates at fermionization resulting in 
the observed notch at the center, in the case of dipolar bosons there is no saturation of the 
interaction energy. 

For the corresponding two-body density in Fig. [2] we observe a gradual depletion of the 
diagonal as d increases from d = and a correlation hole at X\ = X2 can be observed for 
d = 1.4 (Fig. [2] (b)) This is understandable since the particles in the same well attempt to 
minimize the energy by avoiding each other and thus the overlap at x\ = x<i is reduced. For 
stronger interactions (Fig. [2](c)), the diagonal depletion is accompanied by a fragmentation 
of the off-diagonal densities revealing an effective 'fermionization' like pattern. 

B. Incommensurate Filling 

We now study cases with incommensurate filling, which show a much stronger effect of the 
long-range dipolar interactions as well as a stronger dependence on the number of bosons. 
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Figure 1: (color online) One body density for six bosons in the triple well for different interaction 
strengths d. 
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Figure 2: (color online): Two body density for six bosons in the triple well for (a) d = (b) d = 1.4 
(c) d = 15. 

The most important aspect of incommensurate filling compared to that of the commensurate 
case is the absence of a true Mott-Insulator state since there is always a superfluid fraction 
on top of a Mott-Insulator phase. For filling factor v < 1 (two bosons in a triple well) as 
the dipolar interaction increases we observe in Fig. [3] that the atoms are gradually shifted 
to the outer wells. The population of the middle well totally vanishes, in contrast to the 
case of contact interactions [32] • The long-range dipolar forces act here not only locally or 
between next neighbours but also between next to next sites. 

Another strong manifestation of this long-range strong dipolar interaction occurs in the 
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case of a pair of particles on top of a commensurate background (5 bosons in a triple well). 
For intermediate interactions a plateau region in the population (Fig. 4) occurs followed 
by a rearangement of the population with increasing interaction. We observe a general 
upward shift of one-body density at the inter-well space accompanied by a broadening of 
the densities in the outer wells (see d = 1.6 of FigQa)). This is a result of the repulsive 
interaction of the two extra bosons which are pushed to the outer wells. This effect becomes 
striking for higher interaction strengths where the two extra bosons localize in the outer 
wells leading to wiggles in the densities in these wells (d = 2.5 FigQa)). However, as 
the interaction increases further, the intra-well repulsive interaction between the bosons in 
the doubly occupied outer wells forces them to avoid each other, effectively localizing then 
separately in the same well. This leads to the formation of two separate density peaks in 
the one body density of the outer wells as seen for d — 15 (FigQa)). Thus in effect, we 
obtain a 'crystal-like' state consisting of five separate density peaks with two in the outer 
wells and a single peak in the middle well. In a sense the strong dipolar force becomes more 
prominent than the external potential forces, in contrast to the contact interaction which 
acts only locally. 

Fig. [5] presents the two-body densities for different d values. With the transfer of the 
particles to the outer wells due to the long range repulsion and the depletion of the population 
in the middle well, the diagonal contribution is strongly reduced for d = 0.8 (Fig. [51(b) ) . 




With further increase in interaction, we observe the formation of the correlation hole which 
is strongly evident for d = 2.5 (Fig. [5^c)). Moreover, the off-diagonal density peaks starts 
to fragment as the two body interactions gain prominence. We see a complete checkerboard 
formation for d = 15, with depleted diagonal (Figj5^d)) as the strong repulsive interaction 
creates strong fragmentation in the off diagonal densities. 

IV. ENERGY SPECTRUM 

As we shall see later dipolar interactions allows for an interesting multi-band tunneling 
dynamics even for the case of only two bosons. However before we explore the dynamics, 
let us first study the energy spectrum for the two particle case. This will guide us and 
point to the relevant interaction regime and the characteristics of the contributing excited 
eigenstates. 
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Figure 3: (color online) (a) One body density for two bosons in triple well, (b) Variation of 
population in each well with interaction strength d. 

A. Lowest band states 

In Fig. [6j we present the energy spectrum for two bosons as a function of the interaction 
strength d. For weak interactions, the states in the lowest band ( first six states) are 
separated into two groups of triplets. The lowest three are composed of the number-states 
|1, 0, 1), |1, 1, 0) and |0, 1, 1) (the standard number state notation means the population of 
each well eg. the first vector denotes the state where one boson is in the left and one in 
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Figure 4: (color online) (a) One body density for five bosons in a triple well for different interaction 
strengths, (b) Variation of population in each well with varying interaction strength d. 

the right well). However, because of the long-range nature of the dipolar interaction, these 
states are not equivalent energetically ( the |1, 0, 1) state possessing a lower energy) and thus 
the lowest triplet breaks with increasing dipolar interaction. 

The upper triplet is composed of the doubly occupied number-states |2,0,0), |0,0,2) 
and 1 0, 2, 0) with the state |0,2,0) possessing a slightly lower energy because of the hard- 
wall boundary conditions. These states have a much stronger dependence with d as a 
result of the double occupancy in the wells. As the interaction increases, the highest two 
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Figure 5: (color online) Two body densities for five bosons in a triple well for (a) d = 0.0, (b) 
d = 0.8, (c) d = 2.5, (d) d= 15. 

states join together forming a doublet separating from the lower energy state (Fig. [6] inset 
(a)). Unlike the case of contact interactions, here the energy does not saturate in the 
fermionization regime but continues to increase as the interaction strength increases. This 
has as an important consequence for the interplay with the higher band states. 



B. Higher band states 

The increase in energy of the lowest band states with increasing d leads to avoided 
crossings with the energies of the higher band states and to corresponding resonances. The 
higher band states involving primarily singly occupied number-states 1 1 1 , 1° , 0) , 1 1°, l 1 , 0) , 
1 0, l 1 , 1°) , |0, 1°, l 1 ) and 1 1 1 , , 1°) and |1°,0, l 1 ) (where the superscript now refers to the 
band index) have a much weaker dependence of the energy on d compared to the doubly 
occupied number states of the lower band and hence we observe with increasing d avoided 
crossings in the energy spectrum among these states (see Fig. [6] inset (b)). We note that 
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Figure 6: (color online) Two boson energy spectrum as a function of d. Inset: (a) Magnified lowest 
band energy spectrum for small interaction strength, (b) Magnified energy spectrum at the vicinity 
of the avoided crossings. 

since the states |l n , 0, l m ) are energetically detuned from the states |l n , l m , 0) and |0, l n , l m ), 
there are actually two resonances with the lowest band states. These crossings lead to an 
interesting inter-band tunneling dynamics for two particles as we will see next in SecfVj 

In general the two-body spectrum is here much richer than in the case of contact interac- 
tion, since a saturation at fermionization is not forseen for dipolar forces, and thus the long 
range interaction mixes lower and higher states via crossings and anticrossings. 



V. QUANTUM DYNAMICS NEAR RESONANCES 

The avoided crossings occurring in the two-body energy spectrum open the possibility 
for an interesting multi-band quantum dynamics. We thus focus within our study of the 
tunneling dynamics of two bosons on the interaction regime in the vicinity of the avoided 
crossings (d = 0.8 — 1.2). 
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We first consider an initial setup for which both bosons are localized in the left well. 
In Fig. [7] the evolution of the population of the three different wells as a function of time 
is shown for different values of d. For d = 0.8, which is just below the resonance, we 
find a uniform dynamical evolution of the populations possessing a long time period (Fig. 
[7fa)). The initial population in the left well tunnels to the right well with a very long 
period, while the middle well remains almost unpopulated. This is a demonstration of pair- 
tunneling arising primarily from energy resonance of two doubly occupied number-states 
namely |2, 0, 0) and |0, 0, 2). This can be confirmed by studying the pair-probability, that is, 
the probability of finding two atoms in the same well. Fig. [8]^a) shows the evolution of the 
pair-probability and we observe that it remains and oscillates close to unity - implying that 
the bosons tunnel as pairs and there is minimal single particle tunneling. This repulsively 
bound pair phenomenon (or doublon) is also seen and observed experimentally [35J for 
contact interactions, we verify here that it is also a possible mechanism for dipolar forces. 




5000 6000 

Time 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

Time 



Figure 7: (color online) Time evolution of the population density in the left, middle and right well 
for (a) d = 0.8, (b) d = 0.9, (c) d = 0.95, (d) d = 1.0 

For d = 0.9 the resonant region (in the spectrum) is encountered, and we observe a 
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Figure 8: (color online) Time evolution of the pair probability for the same values of d as in Fig. 
7. 

tunneling dynamics containing a combination of many frequencies (Fig. utb)), i.e., many 
eigenstates contribute equally to the dynamics. For long time scales, we observe a dominant 
tunneling from the left to the right well as indicated by their respective populations. However 
now, there is a substantial population in the middle well which oscillates around 20%. This 
long time oscillation is modulated by combinations of fast high frequency oscillations. The 
mechanism of the tunneling could be understood as follows: The tunneling primarily involves 
transfer from |2, 0, 0) to |0, 0, 2) via the intermediate number-states |l n , l m , 0), 0, l m ) and 
|0, l n , l m ) (where m and n indicate the band index). Away from the avoided crossings, the 
intermediate states are energetically far detuned from the initial and final states resulting in 
very long tunneling timescale and minimal contributions of the intermediate number-state. 
Near the avoided crossings, the intermediate singly occupied number state of the excited 
band are energetically in resonance with the initial and final states. This result in both 
faster tunneling and significant contributions of the intermediate states resulting in both 
enhanced single particle tunneling and middle well population. The predominance of single 
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particle tunneling is reflected in the evolution of the pair-probability (Fig. j8j^b)) which 
oscillates between unity and 0.2 implying that in the course of the tunneling it is highly 
probable to find two particles in separate wells, thereby breaking the doublon dynamics. 

Increasing the interaction (d = 0.95) leads to further shortening of the tunneling period 
with even larger modulations due to the participating different frequencies making the pat- 
tern of the tunneling highly irregular as seen in Fig. [7^c). Here, we observe the absence of 
a complete left to right tunneling even for long timescales. Instead, especially for shorter 
timescales, the single particle tunneling dominates. About 2/3 of the population tunnels 
from the left well to the right one with a portion remaining in the middle well, while the rest 
remains in the left well. The mechanism for this tunneling is primarily the resonance between 
the states |2, 0, 0) and |l n , 0, l m ) ± |0, l n , l m ). This is reflected in the pair-probability (Fig. 
[8]^c)) which shows oscillations between unity and near zero, thus displaying the oscillation 
between a pair-state and a singly occupied state. 

A different mechanism is encountered as the interaction is further increased to d = 1.0. 
In Fig. we observe a much larger population in the middle well in the course of the 

dynamics. A complete tunneling from the left to the right well occurs for long timescales. 
Moreover, the numerous large frequency modulations seen in the previous case are here 
reduced although the periodic nature of long-time tunneling is not prominently displayed. 
To understand this we note that unlike that of contact interaction, for dipolar bosons, the 
state |l ra ,0, l m ) becomes energetically offset from the other singly occupied number state 
as a consequence of the long-range nature of the interaction. Hence, a second resonance is 
encountered where the contribution of the state |l n , 0, l m ) becomes minimal. The mechanism 
is primarily a resonance between the |2,0,0) and |0, l n , l m ) states leading to a population 
transfer to both the middle and right well. Also involved in this resonance, we have the 
tunneling channel between the two doubly occupied number-states |2, 0, 0) and |0, 0, 2) which 
is valid especially for longer timescales. The pair-probability (Fig. |8](d)) captures this 
mechanism, showing a two mode oscillation. A very fast oscillation, ranging from unity to 
0.2 implying the breaking of pair and a longer oscillation envelope constitute the dynamics. 
Increasing d further, we go beyond the crossing region thus resulting in a reoccurence of pair 
tunneling (not shown here). 
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VI. CONCLUSION 



We have investigated the ground state properties and quantum dynamics of repulsively 
interacting dipolar boson in a triple well trap. The crossover of the ground state from 
the weak to strong interaction regime has been covered for different filling factors focusing 
primarily on dipolar effects which occur due to the long-range nature of this interaction po- 
tential and differ from that of the pure contact interaction case. In the case of commensurate 
filling, apart from the usual Superfluid-Mott insulator transition, for very strong interaction, 
we observed intra-site localization and extremely strong fragmentation. For incommensu- 
rate filling, additionally we find a general tendency for avoiding the middle well since the 
long range nature of the dipolar force pushes the atoms to the outer wells while the exact 
ground state configuration depends strongly on the exact number of particles considered. 
For strong dipolar forces a crystal like phase appears which is largely unaffected by the 
external potential. 

The two-boson energy spectrum shows avoided crossings between the energy levels of 
different bands. In the region of these crossings, we observe an inter-band tunneling dynamics 
arising from the energetic resonances between the localized number-states corresponding to 
different bands. This regime exhibits different tunneling mechanisms ranging from complete 
pair tunneling to single particle higher band dynamics. This study of the static and dynamic 
properties of ultracold dipolar boson for small bosonic ensembles can serve as a very good 
starting-point for our understanding of the strong interaction regime in larger systems. It 
might also help to develop schemes for creation of novel entangled few-body states and their 
controlled quantum dynamics and transport. 

Acknowledgments 

B. C. acknowledges support from the Landesexzellenzinitiative Hamburg, which is fi- 
nanced by the Science and Research Foundation, Hamburg and supported by the Joachim 
Herz Stiftung. L. C. and P. S. gratefully acknowledge financial support by the Deutsche 
Forschungsgemeinschaft (DFG) in the framework of the collaborative research center SFB 



14 



925. 



[1] M. A. Baranov, Phys. Rep. 464, 71 (2008). 

[2] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. Pfau, Rep. Prog. Phys. 72, 126401 
(2009). 

[3] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, T. Pfau, Phys. Rev. Lett. 94, 160401 (2005). 
[4] Q. Beaufils et al, Phys. Rev. A 77, 061601 (2008). 

[5] M. Lu, N. Q. Burdick, S. H. Youn, B. L. Lev, Phys. Rev. Lett. 107, 190401 (2011). 

[6] K. K. Ni et al, Science 322, 231 (2008). 

[7] S. Yi and L. You, Phys. Rev. A 63, 053607 (2001). 

[8] L. Santos, G. V. Shlyapnikov, P. Zoller and M. Lewenstein et al, Phys. Rev. Lett. 85, 1791 
(2000). 

[9] K. Goral K. Rzazewski, and T. Pfau, Phys. Rev. A 61, 051601 (2000). 
[10] C. Eberlein, S. Giovanazzi, and D. H. J. O' Dell, Phys. Rev. A 71, 033618 (2005). 
[11] K. Goral, M. Gajda, and K. Rzazewski, Phys. Rev. A 66, 023613 (2002). 

[12] T. Koch, T. Lahaye, J. Metz, B. Frohlich, A. Griesmaier and T. Pfau, Nature Phys. 4, 218 
(2008). 

[13] K. Goral, L. Santos, and M. Lewenstein., Phys. Rev. Lett. 88, 170406 (2002). 

[14] C. Menotti C. Trefzger and M. Lewenstein., Phys. Rev. Lett. 98, 235301 (2007). 

[15] B. Capogrosso-Sansone C. Trefzger, M. Lewenstein, P. Zoller, and G. Pupillo, Phys. Rev. Lett. 

104, 125301 (2010). 
[16] G. Bruun and E. Taylor, Phys. Rev. Lett. 101, 245301 (2008). 
[17] N. Cooper and G. Shlyapnikov, Phys. Rev. Lett. 103, 155302 (2009). 

[18] A. S. Arkhipov, G. E. Astrakharchik, A. V. Belikov, and Y. E. Lozovik, JETP Lett. 82, 39 
(2005). 

[19] R. Citro, E. Orignac, S. De Palo, and M. Chiofalo, Phys. Rev. A 75, 51602(R) (2007). 

[20] S. De Palo, E. Orignac, R. Citro, and M. Chiofalo, Phys. Rev. B 77, 212101 (2008). 

[21] P. Pedri, S. De Palo, E. Orignac, R. Citro, and M. Chiofalo, Phys. Rev. A 77, 15601 (2008). 

[22] E. G. Dalla Torre, E. Berg, and E. Altman, Phys. Rev. Lett. 97, 260401 (2006). 

[23] X. Deng and L. Santos, Phys. Rev. B 84, 085138 (2011). 

15 



[24] S. Zollner, G. M. Brunn, C. J. Pethick, and S. M. Reimann, Phys. Rev. Lett. 107, 035301 
(2011). 

[25] M. Maik, P. Buonsante, A. Vezzani, and J. Zakrzewski, Phys. Rev. A 84, 053615 (2011). 

[26] T. Lahaye, T. Pfau and L. Santos, Phys. Rev. Lett. 104 170404 (2010). 

[27] D. Peter, K. Pawlowski, T. Pfau, K. Rzazewski, |arXiv:1201.2615V l (2012) 

[28] B. Xiong, J. Gong, H. Pu, W. Bao and B. Li, Phys. Rev. A 80, 053602 (2009). 

[29] M. Abad, M. Guilleumas, R. Mayol, M. Pi and D. M. Jezek, Phys. Rev. A 81, 043619 (2010). 

[30] H.-D. Meyer, U. Manthe, and L. S. Cederbaum, Chem. Phys. Lett. 165, 73 (1990). 

[31] M. H. Beck, A. Jackie, G. A. Worth, and H.-D. Meyer, Phys. Rep. 324, 1 (2000). 

[32] I. Brouzos, S. Zollner and P. Schmelcher, Phys. Rev. A 81, 053613 (2010) 

[33] L. Cao, I. Brouzos, S. Zollner and P. Schmelcher, New J. Phys. 13, 033032 (2011) 

[34] C. Chin, R. Grimm, P. Julienne and E. Tiesinga, Rev. Mod. Phys. 82, 1225 (2010); 

[35] K. Winkler, G. Thalhammer, F. Lang , R. Grimm, J. H. Denschlag, A. J. Daley, A. Kantian, 
H. P. Biichler and P. Zoller, Nature 441, 853 (2006). 



16 



